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Abstract 

The flux of neutrino-induced muons has been measured with the MACRO 
detector. Different event topologies have been detected, due to neutrino interac- 
tions in the apparatus and in the rock below it. The upward-throughgoing muon 
sample is the larger one and is generated by neutrinos with an average energy 
of ~ 100 GeV. The observed upward-throughgoing muons are 26 % fewer than 
expected and the zenith angle distribution does not fit with the expected one. 
Assuming neutrino oscillations, both measurements suggest maximum mixing 
and Am 2 of a few times 10~ 3 eV 2 . The other event categories due to inter- 
actions of low-energy neutrinos (average energy ~ 4 GeV) have been recently 
studied and the results of these new analyses are presented for the first time 
at this workshop. These data show a regular deficit of observed events in each 
angular bin, as expected assuming neutrino oscillations with maximum mixing, 
in agreement with the analysis of the upward-throughgoing muon sample. 



1 Introduction 

The interest in precise measurements of the flux of neutrinos produced in cosmic ray 
cascades in the atmosphere has been growing over the last years due to the anomaly 

1 The complete collaboration list is given in the paper by G. Battistoni in these Proceedings. 



in the ratio of contained muon neutrino to electron neutrino interactions. The obser- 
vations of Kamiokande (Hirata et al., 1992), 1MB (Casper et al., 1991; Becker-Szendy 
et al., 1992) and Soudan 2 (Allison et al., 1997) are now confirmed by those of Su- 
per Kamiokande (Fukuda et al., 1998) and the anomaly finds explanation in the scenario 
of neutrino oscillation. 

The effects of neutrino oscillation have to appear also in the higher energy ranges, as 
reported by MACRO. The flux of muon neutrinos in the energy region from a few GeV 
up to hundreds of GeV can be inferred from measurements of upgoing muons (Ahlen 
et al., 1995, Ambrosio et al., 1998, 2nd reference). As a consequence of oscillation, the 
flux of upgoing muons should be affected both in the absolute number of events and in 
the shape of the zenith angle distribution, with relatively fewer observed events near 
the vertical than near the horizontal due to the longer pathlength of neutrinos from 
production to observation near the zenith. 

Furthermore the flux of atmospheric muon neutrinos in the region of a few GeV 
can be studied looking at muons produced inside the detector and muons externally 
produced and stopping inside it. If the atmospheric neutrino anomalies are the result of 
neutrino oscillations, it is expected a reduction in the flux of upward-going low-energy 
atmospheric neutrinos of about a factor of two, but without any distortion in the shape 
of the angular distribution. 

Here the measurement about high energy muon neutrino flux is presented, together 
with first results on low-energy neutrino events in MACRO. 

2 The MACRO detector 

The MACRO detector (Ahlen et al., 1993) is located in the Gran Sasso Laboratory, 
with a minimum rock overburden of 2700 hg/cm 2 . It is a large rectangular box (76.6 x 
12 x 9.3 m 3 ) divided longitudinally in 6 supermodules and vertically in a lower and an 
upper part, called attico. The active elements (see Fig. |I|) are liquid scintillator counters 
for time measurement and streamer tubes for tracking, with 27° stereo strip readouts. 
The lower half of the detector is filled with trays of crushed rock absorber alternating 
with streamer tube planes, while the attico is hollow and contains the electronics racks 
and work areas. The intrinsic angular resolution for muons typically ranges from 0.2° 
to 1° depending on the track length. This resolution is lower than the angular spread 
due to multiple scattering of muons in the rock. The scintillator system consists of 
horizontal and vertical layers. Time and position resolution for muons in a scintillator 
box are about 0.5 ns and 11 cm, respectively. 

Thanks to its large area, fine tracking granularity and electronics symmetry with 
respect to upgoing and downgoing flight direction, the MACRO detector is a proper 
tool for the study of upward-travelling muons, generated by external interactions. Its 
mass permits also to collect a statistically significant number of neutrino events due to 
internal interactions. 

3 Neutrino events in MACRO 

Fig. [I] displays the different kinds of neutrino events here analyzed. Most of the detected 
particles are muons generated in Charged Current interactions. Figure |2| shows the 
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Figure 1: Cross view of the detector and topology of the neutrino induced events. The 
stars, the dashed lines and the full lines indicate scintillator hits, neutrino paths and 
charged particle paths, respectively. 



parent neutrino energy distribution for the different event topologies : 

f . Up Through - These tracks come from interactions in the rock below MACRO 
and cross the whole detector (E^ > 1 GeV). The time information provided by 
scintillator counters permits to know the flight direction (time-of-flight method). 
The data have been collected in three periods, with different detector configura- 
tions. In the first two periods (March 1989 - November 1991, December 1992 - 
June 1993) only lower parts of MACRO were working (Ahlen et al., 1995). In 
the last period (April 1994 - November 1997) also the attico was in acquisition. 

2. In Up - These partially contained events come from v interactions inside the 
apparatus. Also in this case the time-of-flight method is applied to identify the 
events, thanks to the attico scintillator layers. Hence only the data collected with 
the attico (live-time ~ 3 years) have been used in this analysis. About 13 % of 
events are estimated to be induced by Neutral Currents or v e CC interactions. 



3. Up Stop + In Down - This sample is composed by two subsamples : exter- 
nal interactions with upward-going track stopping in the detector (Up Stop), 
neutrino-induced downgoing tracks with vertex in lower part of MACRO (In 
Down). These events are identified by means of topological criteria. The lack of 
time information prevents to distinguish the two subsamples. Anyway an almost 
equal number of Up Stop and In Down is expected if neutrinos do not oscil- 
late. Neutral Currents and u e CC interactions constitute ~ 10 % of the sample. 
The analyzed data have been collected with the whole detector with an effective 
live-time of ~ 3 years. 
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Figure 2: Distributions of the parent neutrino energy giving rise to different kinds of 
events, estimated by means of Monte Carlo simulation using the same cuts applied to 
the data. The average energy is ~ 100 GeV for Up Through sample and ~ 4 GeV 
for In Up and Up Stop + In Down samples. 



4 Analysis procedure and results 

The time-of-flight method uses the formula 

1 = cx(T 1 -T 2 ) 
(3 L 

where T\ and T 2 are the times measured in lower and higher scintillator planes, respec- 
tively, and L is the path between the two scintillators. Therefore 1/(3 results roughly 
+1 for downgoing tracks and —1 for upgoing tracks. Several cuts are imposed to re- 
move backgrounds from radioactivity and showering events which may cause failure 
in time reconstruction. Another cut is applied to the Up Through sample requiring 



(1) 



the crossing of at least 200 g/cm? of material in the apparatus in order to reduce the 
background due to low-energy charged upgoing particles produced at large angles by 
downgoing muons (Ambrosio et al., 1998). After all analysis cuts the signal peaks with 
1/(3 ~ —1 are well isolated for the first two samples (see Fig. H). 




Figure 3: 1/(3 distributions after all analysis cuts. In the first plot the Up Through 
data collected in the third period (lower MACRO + attico ) are shown, the shaded 
part indicates events whose 1/(3 value is calculated by means of a fit of three time 
measurements on different scintillator planes. In the second plot the whole In Up 
sample is shown. 

The Up Stop + In Down events are identified via topological constraints. The 
main requirement is the presence of a reconstructed track crossing the bottom scintil- 
lator layer. All the track hits must be at least 1 m far from the supermodule walls. 
The criteria used to verify that the event vertex (or [i stop point) is inside the detector 
are similar to those used for the In Up search. The probability that an atmospheric 
muon produces a background event is negligible. To reject ambiguous and/or wrongly 
tracked events which survived automated analysis cuts, real and simulated events were 
randomly merged and directly scanned with the MACRO Event Display. 

Expected rates and angular distributions have been estimated assuming the atmo- 
spheric v flux calculated by the Bartol group (Gaisser and Stanev, 1995, Agrawal et 
al., 1996). The estimate of v cross-section was based on parton set S\ from (Morfin 
and Tung, 1991), taking into account also low-energy effects (Lipari et al., 1995) for 
In Up and Up Stop + In Down samples. The propagation of muons through the 
rock was taken from (Lohmann et al., 1985). The uncertainty on the expected muon 
flux is estimated 17 % for Up Through events and 25 % for the other events. The ap- 
paratus and the data acquisition are fully reproduced in a GEANT (Brun et al., 1992) 
based Monte Carlo program and the simulated data are processed by means of the 
same analysis chain used for real data. Particular care has been taken to minimize the 
systematic uncertainty in the detector acceptance simulation. For the Up Through 



sample, several different analyses and acceptance calculations, including separate elec- 
tronic and data acquisition systems, have been compared. For each sample, trigger and 
streamer tube efficiency, background subtraction, effects of analysis cuts have been in 
detail studied. The systematic error on the total number of events due to the accep- 
tance has been estimated 6 % for Up Through sample. The uncertainty is higher 
(10 %) for low-energy samples because it depends strongly on data taking conditions, 
analysis algorithm efficiency and mass of the detector. 

In the Up Through sample 479 events are in the signal range (0.25 around 
1/(3 = —1). After the subtraction of the estimated backgrounds, the observed number 
of events becomes 451. For this sample 612 events are expected and the ratio observa- 
tion/expectation is reported in Table [jj Fig. |] shows the zenith angle distribution of 
the measured flux compared with the expectation. The error bars on the data show the 
statistical errors with an extension due to the systematic errors, added in quadrature. 
The observed zenith distribution does not fit well with the expectation, giving a maxi- 
mum x 2 probability of only 0.1 %. The observed number of events and the shape of the 
zenith distribution can be explained in the scenario of — > v T oscillation with a best 
fit point in the unphysical range (sin 2 29 m i X > 1). Both measurements independently 
yield very close mixing parameter values. The second plot in Fig. || shows the best 
fit point in the physical region (sin 2 2# m j Z = 1 and Am 2 = 2.5 x 10 -3 eV 2 ) obtained 
combining event number and angular shape analysis. The probability associated to this 
point is not so high (17 %) because the probability of the zenith distribution is still 
low (5 %), due to the relatively few events in the region — 1 < cos6 < —0.8 compared 
with the number of events in —0.8 < cosO. In the second plot of Fig. ^ the solid lines 
show the contours for 10 % and 1 % of the best fit probability. The dashed lines show 
the allowed regions at 90 % and 99 % c.l. evaluated according to the prescription by 
(Feldman and Cousins, 1998). The dotted line shows the sensitivity which is the 90 % 
c.l. contour which would result from the preceding prescription if the data and the 
M.C. prediction happened to be in perfect agreement at the best fit point. 

The ratios of the observed number of events to the expectation and the angular 
distributions of the In Up and Up Stop + In Down data samples are reported 
in Table |l] and in Fig. [|. The low-energy samples show an uniform deficit of the 
measured number of events over the whole angular distribution with respect to the 
predictions based on the absence of neutrino oscillations. We note a good agreement 
between the results for low-energy and Up Through events. Assuming the oscillation 
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Up Through 


0.74 


0.036 


0.046 


0.13 


1.05 


In Up 


0.59 


0.06 


0.06 


0.15 


1.02 


Up Stop + In Down 


0.75 


0.07 


0.08 


0.19 


0.98 



Table 1: Ratios of observed on expected number of events for different event topologies. 
In the last column the ratio is calculated assuming neutrino oscillation with the param- 
eters suggested by the Up Through sample (sin 2 29 m i X = 1, Am 2 = 2.5 x 10~ 3 eV 2 ). 



parameters suggested by higher energy sample, it is expected a ~ 50 % disappearance 
of Vfi in In Up and Up Stop samples because of the neutrino path (thousands of 
kilometres). No flux reduction is instead expected for In Down events whose neutrino 
path is of the order of tens of kilometres. The ratios and the angular distributions 
estimated assuming the v oscillation are also reported in Table |I] and in Fig. [5| 




Figure 4: In the first plot it is displayed the comparison between measured and expected 
fluxes, for the Up Through sample with a muon energy threshold ofl GeV . The solid 
curve and the shaded region show the expectation for no oscillation and its uncertainty. 
The dashed line shows the predicition assuming v oscillation with maximum mixing 
and Am 2 = 2.5 x 10~ 3 eV 2 . In the second plot the results of the combined analysis for 
Up Through sample are shown in terms of oscillation parameters (see the text). 



5 Conclusions 

The Up Through sample shows a deficit of the measured number of events with 
respect to the predictions based on the Bartol flux in absence of neutrino oscillations. 
The deficit is higher near the vertical direction. Hence the previous results of the 
MACRO experiment (Ahlen et al., 1995) are confirmed. A new paper about this item 
has been submitted for pubblication (Ambrosio et al., 1998, 2nd reference). 

Also the low-energy neutrino events are fewer than expected and the deficit is quite 
uniform over the whole angular range. The three data samples are in agreement with a 
model of —* v T oscillation with maximum mixing and Am 2 of a few times 10 -3 eV 2 . 
The combined analysis of the three different data sets is in progress. 
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Figure 5: Comparison between measured and expected number of low-energy events 
versus cos9 (In Up in the first plot, Up Stop + In Down in the second one). The 
dashed line is obtained assuming neutrino oscillation with the parameters suggested by 
Up Through sample. In the second plot the absolute value of cosO is used because the 
flight direction is unknown. 
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